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We present a differential equation-based mathematical model of nectar foraging by the 



and nectar receivers. Results generated from the model are used to demonstrate how 

different classes within a collective can collaborate to combine information and produce 

finely tuned decisions through simple interactions. In particular we show the importance 

of the 'search time' - the time a returning forager takes to find an available nectar receiver 

- in restricting the forager population to a level consistent with colony-wide needs. 
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1. Introduction 

Collective decisions may result from the simple aggregation of individual preferences. Of- 
ten, however, they emerge as individuals' preferences and their interactions with others 
combine in complex ways. In some important situations, individuals can be categorised 
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by their behaviour into different classes whose interactions drive the decision making pro 



cess. Such scenarios occur in both human and an imal societies ( IConradt and List 



For e xample, social choice theory (jArrow 



19631 ) and its related fields (JDryzek and List 



20091). 



2002) are primarily concerned with how human societies make decisions when classes of 
individuals assign diffe rent utility rn e asures to t he available op t ions. Examples in animal 



societies are given by 



Couzin et al. 



Raghib et al. 



(120101 ). who divide social 



(120051 ) and 

groups such as fish, insects and birds into classes of informed and uninformed individuals 
which each move according to different r ules whilst interacting to ensure the group travels 



cohesively, and by 



Clutton-Brockl ( 11998! ) who discusses how groups decide which members 



will reproduce in societies containing dominant and subordinate females. 



Most collective decision models focus on either the individual-sea 



Pratt et al. 



2005! ) or the group- scale (JNevai et al. 



2009 



Tachikawal . 



e (JBonabeaul . l2002l : 



20101), and therefore 



do not explicitly include interactions between classes which occur at an intermediate level 
of organisation. In this paper we look at the intermediate, or mesoscale, interactions 
between two behavioural classes of Apis mellifera (honeybees), foragers and receivers, 
and explore how this interaction mediates the decision 'how many foragers should we 
send out of the hive to gather nectar?' A representation of these three scales, and how 
they relate to each other, is given in Fig. [H 



Nectar foraging i s a par t icular 



(JLindauer, 



1961 



Seeley 



1995 



y importari t and well studied task of the honey bee colony 



Dyer 



20021 ). Colonies require energy in the form of nectar 
to survive and thrive. Foragers exploit nectar sources (flowers) by collecting their nectar 
and returning it to the hive. The nectar is delivered to waiting receiver bees who then 
process and store the nectar within the hive. The forager's main behavioural choices are 
resting, returning (to the forager source) or re cruiting (performing the waggle dance t o 



encourage other resting foragers gather nectar) ( jCamazine and Sneyd 



1991 



Dyerl . 



2OO2I ). 



Foraging is not only important; it is dangerous and costly. There is a high risk to the 




Figure 1: The level of organisation at which a collective decision is formed may occur 
at different scales. In some cases considering the group as a single entity is appropriate 
(inner circle), in other cases focusing on each of the potentially thousands or millions of 
individuals is more natural (outer circle). For nectar foraging by colonies of A. mellifera 
modelling interaction at the mesoscale (middle circle), that is between a few different 
behavioural classes, allows us to best understand the forces driving the decision making 
process. 



forag er of predation and misadventure resulting in death (JDukasl . 1200 ll : iPage and Peng 



20011 ) ■ A forager expends energy when she flies out on a foraging expedition, which may 
be unsuccessful. Both forager death and unsuccessful foraging may diminish nectar inflow 
and, if frequent, potentially lead to a rapid decline in the colony population (Khoury et al. 
2010) [ submitted - will update upon acceptance or rejection. ]. Clearly, if it is to prosper, 
the colony must make a wise collective decision about how best to allocate foragers for 
maximum nectar return. The principle motivation of this paper is to model the effect 
on forager recruitment of forager and receiver interactions, and thereby understand more 
about their impact on the forager allocation decision. 



Exist i ng models of nectar foraging fo c us, unsurprisingly, on 



1991 



de Vries and Biesmeijer 



1998 



Cox and Myerscough 



f orage rs (ICamazine and Sneydl . 



20031 ). In particular, they fo- 



cus on the recruitment of new foragers through the waggle dance and the communication 
of nectar source quality to other foragers. All such models relate recruitment rates to nec- 
tar source quality; under normal conditions high quality sites will induce extensive waggle 
dancing and high recruitment levels whereas low quality ones may not even lead to waggle 
dancing at all. These models explain how a colony allocates foragers in accordance with 
the quality of the environment, but they do not consider the extent of the colony's need 
for nectar. Through our model we seek to explain how, by balancing quality and need, 
the colony makes a decision that will best ensure its future survival. 



Experienced foragers know about the quality of the foraging environment, which includes 
not only the concentration of sugar in the nectar but also all other factors that a forager 
evaluates when determining the number of waggles runs they will perform, but they know 
little about the state of the hive's nectar reserves. Conversely, receivers know little of the 
quality of the fora ging environraent 



they store nectar (ISeeley et al 



_ Dut learn about the state of the colony's reserves as 
19961 ). We introduce a model that incorporates both for- 
agers (knowledgeable about the external environment) and receivers (knowledgeable about 
the internal environment) to show how interactions between these classes determine the 



size of the active forager population. The model uses a system of three nonlinear ordinary 
differential equations representing the populations of active foragers, foragers unloading 
their nectar, and receivers available to unload nectar. We find both analytic and numerical 
results of the model and use these results to draw conclusions about A. mellifera foraging 
in particular and collective decision making in general. Finally, we extend the model to 
incorporate multiple, dynamic foraging sources and changes in receiver population. These 
models are neither individual-oriented nor group-oriented. They operate at the mesoscale 
(Fig. [1]) in which classes of individuals (in this model foragers and receivers) are the 
entities that are considered. 



2. Simple Model 

We consider two behavioural classes of A. mellifera, foragers and receivers, and construct 
a model to explore how the interaction of these classes determines the active forager popu- 
lation. We assume that the foragers' tendency to rest, to return or to recruit is determined 



by two factors; the qual ity of the forage source 



searching for a receiver (ISeeley and Tovey 



(jSeeley et al. 



199ll ) and the time spent 



19941 ). Initially, we introduce two simplifying 



assumptions; the total receiver population is considered to be constant throughout the 
model and there is a single nectar source of constant quality. Later we will extend the 
model and drop these assumptions. 

The interactions in this model are illustrated in Fig. |2l Two separate sets of data are 
perceived by different classes of the colony, nectar source quality by the forager class and 
colony nectar need by the receiver class. The colony nectar need determines the population 
of receivers and the 'ready to receive' subset of this population. This population, and the 
population of those foragers ready to unload, determines the average time that a forager 
spends searching for a receiver to unload nectar, known as the 'search time'. A forager's 
tendency to choose to rest, to return or to recruit is dependent upon the nectar source 



quality and the search time (ISeeley and Toveyl . [l99J). The outcome of each individual 
decision then feeds back into the population of foragers, and hence into future search 
times. This system places search time at the centre of a web of interactions, and so we 
will describe and define it in more detail. 
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Figure 2: The decision making process of a colony of Apis mellifera engaged in nectar 
foraging. Nectar source quality and colony nectar need are the two inputs, perceived and 
analysed by foragers and receivers respectively. The interactions of foragers and receivers 
produce a value, search time. Nectar source quality and search time are used by foragers 
to decide whether to rest, return or recruit. 



2.1. Search time of returning foragers for receivers 



Upon entering the hive after a successful foraging expedition, a returning forager seeks 
a receiver who will accept her load of nectar. The expected time sp ent searching for a 



receiver can be modelled using the proportion of foragers and receivers ( ISeeley and Toveyl . 



19941 ). Let F be the total active forager population, R be the total receiver population, 
F be the population of foragers seeking to unload near the hive entrance and R be the 
population of receivers ready to receive. Then, assuming that the population in the 



hive entrance is exclusively comprised of foragers seeking to unload and receivers willing 
to receive nectar loads, the probability that a forager's approach to another bee in the 
entrance is to a receiver can be modelled by R/{F + R). If we now let Sg be the constant 
time taken when approaching another bee then the expected search time at time t, S{t), 
is 

Sit) = s,iFit) + R{t))/R{t). (1) 

The expected search time S{t) is critical; it aggregates the outcomes of numerous interac- 
tions between individual foragers and receivers into a single value that strongly influences 
the rate of forager recruitment in the model. It is worth noting that this model, because 
it is a system of ordinary differential equations, can only be used to calculate the average 
search time when, in reality, the actual search time experienced by individual foragers 
will vary significantly around this average. Nevertheless, since we are modelling a colony 
wide phenomenon rather than individual behaviour, using the average search time still 
gives a reasonable indication of the expected behaviour of the system. 



2.2. Equations for forager and receiver numbers 



Let Q be the quality of a nectar source. Here, quality incorporates all criteria that foragers 



assess when determining th eir like 



things like risk of predation ( JDukasl . 



ihood to waggle dance. These criteria might include 



20011) but the most influential component is likely to 
be sugar concentration. Foragers exploiting this source are increasingly likely to recruit 
(i.e. will perform more waggle runs) as Q increases or as S{t) decreases. Conversely, 
their likelihood of resting increases for low Q and for high S{t). Defining fr and fs as 
maximal forager recruitment and resting rates respectively we construct an equation for 
the changes in forager population, F(t), 

dF{t) 



dt 



frF{t) T,{S{t), Q) - f,F{t) T2{S{t), Q). (2) 



recruitment resting 



Here Fi and r2 are functions that scale the recruitment and resting rates by mapping 
source quahty and search time to values between and 1. Fi should increase recruitment 
rates as quality improves or as search time decreases. F2 should increase resting rates as 
quality decreases or as search time increases. We have chosen to define these functions 
using Hill equations since this type of equation permits a good phenomenological fit and is 
algebraically convenient. We define mq and ms as the values that produce half- maximal 
scalings for Q and S{t), and j and k as the coefficients that control the sensitivity of the 
functions to changes in Q and S{t). This leads to definitions for Fi and F2, 



USit),Q)= ^''' ,^ , (3) 



nig Q^ 

m% QJ -] 






In Fig. [3a] we show how different search times affect the search-related part of Fi for a 
range of coefficient values. Compared to foragers, there is a paucity of information about 
rec eivers i n the existi ng literature and so we rely extensively on observations re corded 



bX 



SeelevI ( 1992 



I995I ). We find that the closest comparison to the data given in 



Seeley 



(119921 . Fig. 7) is with k = 4 and ms = 10. Similarly, we plot the quality term of Fi for a 



range o: 



coef ficients in Fig. |3b]aiid determine that the best fit to the data in Fig. 5.31 of 
Seeleyl ( ll995l ) is j = 4 and mq = 1.5. In the absence of similar data for resting rates we 
used these same coefficients for the equivalent terms in F2. This also keeps the number 
of parameters in the model to a minimum. Eq. [2] emphasises nectar source quality and 
search time, it does not consider other factors such as the tendency of foragers to resume 
foraging at sources they were exploiting the previous day. 

Next we define the changes in the population of foragers that are seeking to unload nectar 
within the hive. These foragers, denoted F, are a subset of the total foragers F. The 
number of foragers not seeking to unload, i.e. those outside the hive gathering nectar, is 
F{t) — F{t). These foragers arrive at the hive with nectar at a rate determined by the 
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Figure 3: Plots of the Hill equation terms in Fi with different coef ficient values . Hal f 



maximal values are derived by comparison with experimental data (jSeeleyl . Il992l . Il995l ). 
Coefficient values are chosen by comparing these curves to this experimental data. Fig. 
3al shows how increases in search time decrease recruitment. For very low search times, 
recruitment (effectively the number of waggle runs the forager performs due to quality) is 
essentially unchanged. For very high search times waggle dancing may effectively cease. 
Three different curves are shown representing different values of the coefficient k. Fig. [3b 
shows how increases in quality lead to foragers recruiting more. For low quality sources 
recruitment is negligible. As the quality approaches 1.5, recruitment is half the maximal 
rate. As quality approaches 3 the rate of recruitment plateaus. Three different curves are 
shown representing different values of the coefficient j . 



average foraging trip time fa- They then take an average of S{t) seconds to find a receiver 
and commence unloading. The equation for this population is 

dFjt) _ Fjt) - Fjt) _ m 

dt fa Sit) • ^ ^ 

rate of arriving in entrance rate of finding a receiver 

For this simple model we assume that -R(t), the total number of available receivers, does 
not change and so R{t) = Rq, the initial receiver population. Each receiver can be either 
waiting to receive nectar or engaged in nectar storage tasks. The population of receivers 
engaged in storing is given by R{t) — R{t). They complete this task at a rate determined 
by the average time for nectar storage, denoted r^. Meanwhile, as foragers find receivers 
available for unloading, these receivers drop out of the available receivers population and 
join those storing nectar. The equation for the population of available receivers is 



dR{t) _ Rq- R{t) _ F{t) 

dt ~ Vs ~ Sit)' 

rate of returning from storing rate of receiving nectar 

The complete system of ODEs from Eqs El [3l HJ |5] and [6] is 



(6) 



dFit)_^^.^ m| Q^ ^^,. Sit) 



frFit)-—r^^—^— - fsFit)- 



^ mi 



dt ''' ' ' 5(t)'= + m| QJ + m^ ' ' ' ' Sit)^ + m| QJ + m^ 
dPit) _F(t)-F(t) Fit) 

~^r~ Ja W) 

dRit) _Ro-Rit) Fit) 
dt ~ Vs Sit)' 



(7) 



Unless otherwise stated, we will use the parameters given in Table [H In many cases a 
single value has been chosen from a range of appropriate values (e.g. the average forager 
trip time varies greatly). 
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2.3. Steady states 



We seek information about the structure and stability of the system and so consider 
the equihbrium points {F*,F*,R*) of the simple model (Eq. [Tj). This system has three 
equilibrium points. The first corresponds to the state observed at the start of the day in 
which no foragers are active, the second to the population mix that the colony approaches 
over the course of the day in a normal foraging environment and the third to a biologically 
impossible situation where the population of foragers seeking to unload is negative. Here 
we will restrict our attention to the first two equilibria, the one near which the colony 
typically begins and the one which it typically approaches during the course of a day. 

The first equilibrium point is 



(fA 




( \ 

\ 


F* 


= 





\r-] 




\n.) 



in which there are no foragers active and no receivers storing nectar. To find when the 
forager population increases away from this point we determine when it is unstable. We 
find that this steady state is unstable whenever Q > mQ{fs/ frY^^ (ss/ms)''^^ . For the 
parameter values of Table[T]this corresponds to Q > 0.5, so for all except very low quality 
environments the forager populations will increase from zero. 



Parameter 


Symbol 


Value 






Source 


Single interaction time 


Ss 


5 s 






Seelev et al. fl991") 






(range is 2 s — >■ 


7 


s) 




Forager recruitment rate 


fr 


0.0010 s-i 






Seelev and Towne fl992) 


Forager resting rate 


fs 


0.0002 s-i 






Camazine and Snevd fl991") 


Half-maximal search time 


ms 


10 s 






Scclcv (1992) 


Search time coefficient 


k 


4 






Scclov (1992) 


Half-maximal forage quality 


'"Q 


1.5 






Scclcv (1986) 


Forage quality coefficient 


J 


4 






Seele^ (1986) 


Forager round-trip time 


fa 


15 min 






Seelgi (19941 






(range is 1/2 'min 


-5> 20mm) 




Receiver storage time 


Ts 


20 min 






Seelev (1995) 






(range is 1 m,in 


-s 


20mm) 





Table 1: Standard parameter values 
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The second equilibrium point is 



(fA 


/ 


p* 


= Ro 


\r* ) 


\ 



77 (/a + "^57) 

^7-1 



Ssms-y 



irs+Ss)ms'y-rsS^ 



{rs+Ss)ms-y-rsSs 



\ 



J 



(f-] 




( U \ 




Ro 




F* 


Ts 


ms-f 


\r* j 




[ Ss ) 



for 7 = d yQ^^^~q ■ This shows that the final populations are linearly related to the total 
receiver population, and are only slightly dependent upon source quality (through 7). In 
other words, the population of receivers determines how many foragers should be allocated. 
Calculating the stability of this point is analytically intractable, but numerical simulations 
indicate it is an attracting stable point and an analytic analysis of an approximation to 
it, i.e. of 



(8) 



also suggests it is always stable. This analysis of the system dynamics indicates that in 
a normal foraging environment the forager population will increase from zero towards a 
single stable value. It is surprising that this value is principally dependent on forager 
arrival rate, receiver storage time and receiver population. In some cases, the stable 
equilibrium point in the model is approached but not reached within the limited time 
that a colony has to forage each day. We note that for low search times, ^F{t) ^ 
-^— — '-j-^F{t), so that F{t) Ri FQexp{^—- — ^-j-^t), where Fq is the initial population. 
Clearly, the population grows faster for higher values of Q than for low values of Q. 
So Q strongly determines the rate that the colony approaches the steady state (in the 
model) but only weakly affects the value of that steady state. This is illustrated for the 
parameters of Table [1] in Fig. |H This suggests that the forager population is more likely 
to attain (and even, temporarily, exceed) the stable equilibrium point in a high quality 
environment than in a low quality one. To understand the details of how this works in 
practice, we solve the system numerically. 
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Figure 4: The population of foragers after 4, 8 and 12 hours for simulations with different 
values of Q and otherwise identical parameters. Actual populations are shown as stars, 
squares and diamonds, the equilibrium population as a circle. The equilibrium point 
changes only slightly as quality increases as it is principally dependent upon R. For low 
values of Q (i.e. Q < 1), the population does not grow close to equilibrium, even after 
twelve hours. For moderate values of Q (i.e. Q G [1,2]), the population approaches 
equilibrium quickly and smoothly. For high values of Q (i.e. Q > 2) the population 
rapidly exceeds the equilibrium population and then, upon entering a stabilising phase 
dominated by inhibition of recruitment, slowly decreases towards it. 
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3. Numerical Results 

Fig. [5] displays the results of numerical simulations of Eq. [7] with the parameters of Table 
[H The simulated period was eight hours. This period was chosen because it conforms to 



the periods used in existing models (jCamazine and Sneyd 



1991 



1998 : 



Cox and Myerscough 



de Vries and Biesmeijerl . 



20031 ). Each of the four plots (Fig. [5al Fig. [50 Fig. [5c] and 
Fig. [5d|) show the population timecourses for simulations using a different combination 
of quality and receiver numbers. Fig. [5a] and Fig. [5b] are generated from simulations 
that use the same high quality environment {Q = 3) and vary only in receiver population. 
This gives the simulation results in Fig. [5bl that has half the receiver population of the 
simulations displayed in Fig. [Sa] producing half the forager population. The shape of the 
population timecourse is preserved and the search time is identical. The simulations used 
to produce Fig. [5c] and Fig. [5d]also differ only in receiver populations, however they both 
have a low quality environment {Q = 0.9). Very different results are produced to those in 
Fig. [5a] and Fig. [Sb] the populations are never near equilibrium and, most significantly, 
there are no differences between their population timecourses. In other words, the size of 
the receiver population does not affect the forager population in a low quality environment 
but strongly affects it in a high quality environment. 

These results suggest that, in a high quality foraging environment, the population of 



foragers rapidly increases, as discussed in Sec. 12. 3[ until the number of foragers seeking 

to unload becomes large relative to the number of receivers available for unloading. Once 

this occurs the search time increases and recruitment is inhibited, leading to much slower 

changes in population in accordance with Eq. [2] The transition between these 'rapid 

growth' and 'stabilising' periods can be quite sudden. In Fig. [5a] it can be observed to 

occur at a cusp in the forager population timecourse just before t = 2. This cusp typically 

occurs near or above the equilibrium population of F. In a low quality environment the 

rate of increase in the forager population is far slower and so the the number of foragers 

remains small and search time never becomes large. 
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Figure 5: Timecourses using the model defined by Eq. [7] with varying Q and R and 
other parameters as defined in Table. [TJ Figs. [5a] and [5b] have high quality {Q = 3.0) 
environments and Figs. [5c] and JSd] have low quality {Q = 0.9) environments. Figs. [5a] 
and [Sc] have a large number of receivers {R = 100) and Figs. [5b] and [5d] have a smaller 
number of receivers {R = 50). The legend is defined in Fig. [Sa] Changing R dramatically 
affects the results of simulations using high values of Q, however this has no effect on 
simulations using low values of Q. 



15 



This suggests that in a low quahty environment the population of foragers is highly 
sensitive to Q, which controls its rate of growth. In a high quality environment the rate 
of growth is so fast that the onset of high search times is most important in determining 
the forager population. We find then that the quality, perceived by foragers, determines 
the population in low quality environments. When the quality is high however, it is 
receiver numbers that determine the forager population, since by Eq. [T]the search time is 
dependent on the receiver population. That is, we can consider the quality of the foraging 
environment as providing positive feedback to foragers and the receiver population as 
providing negative feedback to foragers. The positive feedback dominates when the forager 
population is small and the negative feedback dominates when the forager population is 
large. 

We offer a biological interpretation for this. A colony with large nectar reserves may 
need to risk only a small number of foragers in source exploitation. An efficient and 
simple way to do this is to limit the population of receivers, which will increase search 
time once this small number of foragers is active, preventing further recruitment. When 
the colony has depleted its reserves, for example after an extended period of rain, it can 
make many receivers available and foraging numbers will become constrained only by 
forager population size. Foragers, however, need not respond to high receiver numbers 
by recruiting to the maximum level if the nectar source quality is low. Thus, the two 
behavioural classes are interpreting different sets of data, i.e. foraging environment quality 
or colony nectar need, and using these different sets of information to collaboratively reach 
a decision about the appropriate population of foragers. The collaboration is via the search 
time, a value constructed in a self-organised manner through the interactions of foragers 
and receivers. This interpretation is biologically sensible; receivers, who store nectar, 
presumably understand the need of the hive for nectar and hence, by their accessibility 
for unloading, communicate to foragers that the hive does or does not require a greater 
intake. Likewise, foragers have knowledge of the foraging environment and are best suited 
to determining that the quality is too poor to be worth investing in more expeditions. 
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A graphical representation of this is given in Fig. |6l This figure shows the population 
of foragers eight hours after the commencement of foraging against a range of values for 
R and Q. It is clear that, for most values, the population is sensitive to either R or 
Q but not to both. That is, for low qualities it is the quality of the forage source that 
determines the forager population and for high qualities it is the receiver population. For 
Q < 1, the forager population is typically still growing after eight hours. For some values 
of R there is a ridge at Q ~ 1 because the population in the simulation with a slightly 
higher quality has already commenced its stabilising phase and its forager population is 
already decreasing (c.f. Fig. [5a|) . For Q G [1, 2] the plot is fiat along the Q axis since the 
simulations have reached equilibrium within eight hours. For higher qualities the plot's 
steepness increases, refiecting those simulations whose forager populations have peaked 
well above equilibrium and not reached equilibrium after eight hours (c.f. Fig. H]). 

4. Extending the model with receiver recruitment via the tremble dance 



The tremble dance (ISeeley et al. 



1996 



Dyer 



2OO2I ) is performed by foragers experiencing 



high search times, and is associated with increases in the receiver population. We seek to 
understand how receiver recruitment affects the collective decision about forager numbers, 
and so we extend the simple model (Eq. [7]) to include tremble dancing. Let U be the 
unit step function and u, m and rriT be, respectively, the maximal recruitment rate due 
to tremble dancing, the 'steepness' of the response to tremble dancing and the minimum 
search time for which tremble dancing occurs. The equation that we use to model receiver 
population change is 

--li = U(S(t) - my „, , ^ . 9) 

Here there is no change in receiver population unless the search time is greater than rriT- 

When S{t) > ttit the receiver population increases at an increasing rate as search time 

increases. For w = the system reduces to the simple model of Eq. [71 The values of itit 

and m were determined in a similar manner to m^, mQ,j and k, by plotting for a range of 
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Figure 6: Forager populations after eight hours against a range of hkely Q and R values. 
At a threshold of Q ~ 1 the population switches from being sensitive to Q to being 
sensitive to R. The populations may exceed the stable equilibrium point temporarily. 
This is apparent around Q = 1 where the graph displays a dip as Q increases since 
simulations with quality Q ^ 1 have decreased further from their peak after eight hours 
than those with a lower quality. This corresponds to the cusp that occurs as the forager 
numbers change in the transition from the 'rapid growth' to the 'stabilising' phase. A 
similar effect is visible for Q > 2, where there is considerable over-shoot beyond the 
equilibrium population and the populations are yet to decrease to equilibrium. 
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values (Fig. [7]) and comparing to (ISeeleyl . Il992l . Fig. 7); we find that m = 5 and m^ = 30 
are suitable values. We begin our analysis of the new system by calculating the effect that 




50 100 

S(t) 



150 



Figure 7: The effect of different values of m on the trembling response. The best fit ting 
median and coefficient values were chosen by comparing with results in ISeeleyl (119921 ). 



the introduction of Eq. [9] has on the steady state. F{t) is unchanging whenever F{t) = 
(trivially, though this may be unstable when the population is perturbed) or whenever 



S{t) = S* 



IQ^ 



fsm^o 



(10) 



When S* < ijit, Rit) eventually stabilises and then the rest of the system also stabilises. 
If, however, S* > itit then the search time steady state S* (Eq. [TOj) continually induces 
receiver population increases and the entire system loses stability. By rearranging Eq. [10] 
we can see that the system will be stable whenever Q < tuq {i'm!^/mg){fs/ fr)) '' ■ For 
the standard parameters used in the simulations and defined in Table [1] this implies that 
fo^^ Q ^ 3 the forager population should reach a stable value, for higher values of Q it will 
continually increase. 



This stability criterion is independent of -R(t), so although high search times will increase 

the receiver population, this increase does not resolve the instability. Hence for high 

quality sources, the increases in both receivers and forager populations will continue until 

the colony has exhausted its populations of potential foragers or potential receivers. 
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This is observed in the numerical simulations presented in Fig. [HI Fig. |8a] is generated 
using an identical configuration to Fig. [SaJ with the extra inclusion of Eq. |9]with w = 5. 
The forager and receiver populations become unstable, continually increasing. In reality 
the colony would soon exhaust its population of potential foragers. A series of jagged 
oscillations occur for the search time when it exceeds the threshold ttit defined in Eq. 
[TUl When the search time exceeds this threshold there is a sudden infiux of receivers 
which then causes the forager population to increase until the search time threshold is 
reached again. In Fig. |8b] we have a similar set of parameters, however we have reduced 
the foraging quality so that Q = 2.0, which is below the stability threshold. In this 
simulation, as predicted, the forager population trends towards a single value. During the 
rapid growth phase, receiver recruitment occurs due to the high amplitude of the search 
time oscillations, but this ceases as the system stabilises. A simulation (not shown) 
generated over the same parameter set except with Q = 0.9 produces results identical 
to that in Fig. [5cl We find here that the inclusion of Eq. [9] has no effect in this 
low quality environment, since S{t) is always less than m-r- This is consistent with the 
previous conclusions concerning low quality environments in that receiver numbers have 
no influence over the forager population. 

These results describe a mechanism by which the equilibrium population of foragers 
may increase. Due to, amongst other factors, the high degree of genetic diversity in A. 
mellifera colonies, the many stim uli thresholds at which individuals initiate action vary 



(JMyerscough and Oldroyd 



20041 ): this is no doubt true for the tremble dancing response 



amongst the pool of potential receive r bee s. The effect of a mechanism similar to those 



described by 



Fuchs and Moritzl f ll999l ) and 



Jones et al. 



(12004[ 1 would be beneflcial as the 



stable state populations would increase smoothly as the source quality increases. When 
there is a range of response thresholds, receivers play a similar role to that described in 
this simple model. Fundamentally, however, this model states that receivers can adapt to 
exceptionally high quality environments in such a way that forager numbers can increase 
further and thereby take full advantage of foraging conditions. 
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Figure 8: Population tiinecourses for a model with parameters identical to Fig. I5al except 
that Eq. |9] is included in the model with w = 5. The dynamics are identical up to 
t ?« 2, at which point tremble dancing commences and receiver numbers begin to increase. 
This point coincides with the end of the rapid growth phase and commencement of the 
stabilising phase, discussed in Sec. O In Fig. [SalQ = 3 and S* > rriT, and so the receiver 
population increases indefinitely. The search time undergoes oscillations due to repeated 
cycles of receiver increase followed by forager increase causing the populations to vary 
widely. In Fig. |8b] Q = 2 and S* < uit, and so the receiver and forager populations 
eventually stabilise. 
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5. Extending the model with multiple forage sources 

In the field an A. mellifera colony will simultaneously exploit multiple sources, each with 
different quality. We therefore extend the model defined by Eq. [7] to include more than 
one nectar source. We denote the quality of n sources as Qi, Q2, ■ ■ ■ , Qi, ■ ■ ■ , Qn, and add 
new forager equations for each source so that Eq. [7] is replaced by 

^^^^*) ^/,F,(t) T,iS{t),Q,)-fMt) T,iSit),Q,), 



dt 

dF^it) _ F,{t)-F,{t) F,{t) 

dt fa S{t) ' 

dR{t) _R{t) - R{t) EliHt) 



(11) 



dt Ts S{t) 

and search time (Eq. ^ is given by S{t) = s, (^"^^ (^i(t)) + i?(t)) /R{t). 

We use this extended model to investigate how the inhibitory effect of high search times 
affects the population and the allocation of foragers among different sources, and to 
explore whether the results of the single source model (Eq. [7]) also apply to multiple 
source models. 

The results of simulations with two nectar sources, with qualities denoted Qi and Q2, are 

presented in Fig. |9l Fig. [9a] displays the population timecourses from a simulation of a 

high quality environment, in which Qi = 3 and Q2 = 2. The total population of foragers 

in this plot is approximately equal to the population of foragers in Fig. |5al which also 

presents results from a simulation of a high quality, though single source, environment 

with Q = 3. The equilibrium populations of the two simulations are identical. The 

similarity between the total populations is consistent with the idea, explored previously, 

that the receiver population imposes an upper limit on the forager population. The 

simulation suggests that receivers achieve this by limiting the total population of foragers 

rather than by limiting the population at each nectar source. Further, the time when 

the population djTiamics move from a 'rapid growth' phase to a 'stabilising' phase is the 

22 



same. Other characteristics of the figures are very similar, such as search times, available 
receiver population and the shape of the timecourses. For a variety of measures then, the 
single source quality model of Eq. [7]is equivalent to aggregating the population of foragers 
at multiple sources in high quality environments. This is expected because receivers don't 
know which nectar source foragers arrive from. 

The distribution of foragers between the two sources changes during the rapid growth 
phase, with the higher quality source acquiring a larger percentage of the foragers until 
the stabilising phase is entered, when the distribution becomes nearly static. The model 
assumes that there is no mechanism by which receivers direct the foragers to one particular 
source - they do not 'manage' foraging. Instead, they specify the colony-wide nectar goal 
via the receiver population size and foragers determine their allocation to each source 
based on each source's quality. Te model predicts that, as usual, bees will display a 
decentralised and self-organised approach to colony management. 



Fig. [9b] presents the population timecourses from a simulation of a low quality environ- 
ment, in which Qi = 0.9 and Q2 = 0.91. The forager population at each source is equal to 
that of the population of foragers in Fig. [ScJ which presents results from the correspond- 
ing simulation of a single source environment with Q = 0.9, and hence the total forager 
population in this plot is extremely close to twice that of the single source simulation. 
This is consistent with the results of the model with a single low quality source, in that 
foragers determine their own population on the basis of source quality without being in- 
fiuenced by receivers. The distribution of foragers is determined by the quality, and since 
Qi and Q2 are nearly identical so are the populations foraging at these sources. 

The numerical results of the multiple source simulations are consistent with those of the 
single source simulations. This suggests that the analytic results for the simple model 
(Sec. E]) are also relevant to this multiple source model. As such, we can consider the 
single source model as aggregating multiple sources into one general term by thinking 
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of Q as representative of the quality of the entire environment rather than just a single 
nectar source. 



6. Extending the model with dynamic forage sources 



The quality of a nectar source may change over the course of a day (jGinsberg) . Il983l ). 
due to changes in weather, inter-species competition levels or plant characteristics. How 
a colony adapts to this is crucial to its success, and so we seek to understand how search 
time affects f orage r populations when source quality changes. We follow the scenario of 



Seeley et al. 



(119911 ) in which there are two nectar sources, one initially of low quality and 



one initially of high quality, that switch qualities after a few hours (denoted Tg) so that. 



2.5 if t<T, 0.9 ift<T, 

Qi= { and Q2= { 

0.9 otherwise 2.5 otherwise 



Using the simple model of Eq. [TJ we generate simulations of two different scenarios 
with this type of nectar quality switching. Each simulation initially has one hundred 
receivers with other parameters as defined in Table [H Population timecourses of the 
forager populations and search times are presented in Fig. [TOl In Fig. llOal T^ = 1 and 
the foraging popula tion rapidly changes its efforts to exploiting the better source, as with 



Seeley et al. 



(Il99l[ ). In Fig. llObI the parameters are identical but T, = 2. The better 



source is not exploited in this scenario, suggesting a failure in the decision making process. 

The colony's inability to adapt to the quality change in Fig. llObI is surprising given its 
successful adaption to the apparently similar scenario recorded in Fig. llOal The reason 
for the failure is that during the first simulation the switch occurred during the 'rapid 
growth' phase, during which search times are small, and the second simulation's switch 
occurred during the 'stabilising' phase, during which search times are high. High search 

times inhibit recruitment, regardless of source quality, and so during the second simulation 
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Figure 9: Population timecourses for a model with parameters identical to those in Fig. 
[5], except that there are two forage sources. In both plots there are one hundred receivers, 
so they correspond to Fig. [5a] and Fig. [5cl Fig. |9a] has two high quality nectar source 
(Qi = 3,(52 = 2) and Fig. |9b]has two low quality ones (Qi = 0.9, Q2 = 0.91). 
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there simply wasn't enough recruitment occurring to increase the forager population at 
source two. By contrast, in the first simulation there were low search times for more 
than two hours after the source change. Previous results consistently show that the rapid 
growth phase is associated with fast changes in population whereas the stabilising phase 
is associated with slow population changes. We note, further, that by Eq. [21 if 

Fiji) ^ Qi Qi + mj^ 

and Qi < Q2, then ^ > ^ and hence there is still more recruitment to the poor quality 
site than to the high quality one. For the parameter values used in the simulation, this 
means that for ^jn > 15 (i.e. when there are 15 times more foragers to source one than 
to source two) the first source will continue recruiting more foragers even after it becomes 
the poorer quality source. This can be observed in the timecourses of the populations 
in Fig. llObI around t = 4. Nevertheless, the better quality site will increase its forager 
population at a greater proportional rate and so eventually, over a period of time much 
longer than a day, have the greater population as well. 

With such a poor result observed in Fig. I10b| it is natural to ask if there is something 
missing from the model that could improve the outcome and make it more realistic. We 
will apply the same values as used to generate Fig. llOal and Fig. llObI into a model 
extended to include tremble dancing (Eq. |9]) with u = 0.1. The resulting population 
timecourses are displayed in Fig. llOcI and Fig. llOdI The results of Fig. llOcI are similar to 
those in Fig. llOal though with the increased populations expected due to the inclusion 
of tremble dancing. The results of Fig. I10d| however, improve considerably upon those 
in Fig. llObI The improvement occurs because tremble dancing relieves some inhibitory 
pressure on forager recruitment, as the increased receiver numbers temporarily decrease 
the search time. Increasing numbers of these new recruited foragers exploit the higher 
quality nectar source. 
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Even with the modifications, the change in forager allocation when the quality switch 
occurs at Tg = 2 still appears to be slow. There are two possible reasons for this. Firstly, 
the model could be showing a result that is qualitatively accurate, yet its outcome is 
much better in practice due to behaviours that are not modelled. Secondly, the model 
could be quantitatively correct and this slow response is advantageous in someway to the 
colony, representing perhaps a collective memory effect that highly values nectar sources 
that have previously proven to be of good quality. It could also be seen as a part of 
the colony's risk-reduction strategy, as the colony allocates its foragers sub-optimally (in 
terms of maximising nectar influx) so that it can robustly protect against another sudden 
loss of the better source. 



7. Discussion 

We have presented a class-oriented model of A. mellifera nectar foraging whose principal 
novelty is that it incorporates the receiver population. Numerical and analytic results 
derived from the model all suggest that the size of the receiver population determines the 
equilibrium population of active foragers. In general, foragers exploiting a high quality 
nectar source rapidly recruit more new foragers, but only while there are receivers available 
to quickly unload them. Once the demand for receivers significantly exceeds supply, 
forager recruitment is inhibited and their population stabilises. 

Under the simple model defined and analysed in Sec. |2]and Sec. [3l receiver numbers are 

fixed. Once this fixed population is reached no more bees will become receivers (though 

until it is reached some of the receivers may be engaged in other worker tasks). The model 

is modified in Sec. |4]so that new receivers are recruited by tremble dancing. The results 

of the simulations produced by a model that includes tremble dancing suggest that when 

an extremely high quality nectar source is found, the colony may indefinitely recruit new 

receivers, at least until the worker population is exhausted. Normally however, tremble 
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Figure 10: Population timecourses for simulations of switching sources. Sources switch 
quality between 2.5 and 0.9. Parameter values are from Table [1] and R = 100. In Fig. 
llOal nectar sources switch quality at t = 1 and there is no tremble dancing. The colony 
adapts to the changed environment. In Fig. llObI nectar sources switch quality at t = 2 
and there is no tremble dancing. The colony takes a long time to adapt to the changed 
environment. In Fig. llOcI nectar sources switch quality at t = 1 and there is moderate 
tremble dancing. The colony rapidly adapts to the changed environment. In Fig. llOdI 
nectar sources switch quality at t = 2 and there is moderate tremble dancing. The rate 
and extent of adaption is improved as comj^ed with Fig. llObi 



dancing allows the colony to increase its population of receivers to a greater, but still 
limited, size if this is justified by the discovery of a high quahty nectar source. 

Another system is also plausible; some worker bees are dedicated receivers, others are 
available to receive under certain criteria and others are never available. Such a system 
allows a sta ggered response to tremble dancing, so that individual bees may have different 



thresholds (jjones et al 



2004J ) at which they will switch to receiving; some will switch as 
soon as they perceive a tremble dance, others only when the colony need is great. The 
actual scheme is of great importance in determining the role of receivers. If workers 
always switch roles to become receivers upon noticing a tremble dance, then they are only 
passive participants in the decision making process. If, however, they are influenced by 
their environment as well as the tremble dance then they have an important role to play 
in deciding how many foragers should be active (as we have modelled here). The effect 
on the decision making process is immense, and warrants further study. Under such a 
scheme, the tremble dance recruits non-receivers to act as receivers but this will only occur 
amongst those bees who perceive a high colony need. Once again, this makes biological 
sense; bees that are available to receive should be engaged in other housekeeping roles 
until such time as there is demand for their labour. 

The results reported in Sec [5] suggest that receivers are not directly involved in allocat- 
ing foragers to specific forage sources. Instead, their population imphcitly specifies a 
target forager population and the foragers determine the distribution amongst different 
sources. We note, however, that there is experimental evidence that receivers can com- 



pare the quality of nectar coming in from dif ferent sources and pr e 



forag ers carrying the highest quality nectar (jScheiner et al. 



2004; 



erentially respond to 



Martinez and Farinal . 



20081 ). These foragers will have a shorter search time (and a shorter unloading time). 
Receivers may influence the allocation with this mechanism, since recruitment to high 
quality sources will be increased and recruitment to low quality sources decreased as 
foragers returning from high quality sources have a longer search time. 
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Sec. [6] models a colony adapting to a dramatic change in nectar source quality. Interpreted 
qualitatively, the colony response depends on whether the change in quality occurs before 
or after the start of the 'stabilising' phase - i.e. when the demand for receivers outstrips the 
supply. The colony smoothly adapts to a change before this phase, but takes much longer 
to respond correctly when the change occurs afterwards, when recruitment is inhibited by 
high search times. The inclusion of tremble dancing improves the response significantly. 
This is interesting because it highlights how flexibility within the colony, in this case in 
the receiver population, leads to improved outcomes. 

7.1. Experimental implications 

Nectar receivers have received less attention than foragers. However, our results indicate 
that the nectar receiver population has a significant influence on the colony's exploitation 
of forage sources. For instance, the model suggests that the daily foraging cycle may 
consist of periods in which the nectar receiver population inhibits forager population 
change - even as conditions themselves change. It would be interesting to see whether 
this inhibition does in fact occur in real hives. 

In addition, the population of active nectar receivers may fluctuate due to colony need 
and the level of tremble dancing. The model results suggest that receivers should respond 
to tremble dancing only when the foraging environment is excellent and/or the colony 
has low nectar stores. Observations are needed to demonstrate that the activity of the 
nectar receivers indicates the internal state of the colony and that this, in turn, affects 
the forager population. 

7.2. Model limitations and future extensions 

A number of extensions are suggested by the biology and by the model results. Over long 

periods a hive may fill up with nectar. Until new honeycomb is constructed this has the 
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effect of increasing storage time; storage time is represented in this model by the constant 
parameter r^. We see from Eq. [8] that as r^ increases, the steady state populations 
decrease. The model could be extended to make r^ vary with the level of nectar stored 
within the hive. This may result in greater nectar stores leading to smaller foraging levels, 
and exploring this in the model may elucidate the system further. In addition, including 
the tendency of workers to construct new honeycomb as the hive fills and being able to 
visualise its effect on the forager and receiver populations would be useful and interesting. 

Another extension is to allow receivers to respond more favourably to nectar gathered 
from particular locations. We have discussed previously how this could lead to increased 
foraging populations at the favoured nectar sources, but modelling it would quantify this 
and allow us to better understand its potential effect under different circumstances. 

One further extension would be to include inspector bees (foragers who repeatedly check 
a previously good quality nectar source to ascertain whether it has returned to being good 
quality) and different subclasses of foragers that respond differently to different values of 
Q and S{t). This should be particularly relevant when studyi ng task allocation, models 



of which often rely on different individual response thresholds (JBeshers and Fewell 



Myerscough and Oldroydl . 



2001 



2004 . 



7.3. Foraging as a collective decision 

Our model treats foraging as a problem solved by a collective group intelligence; the 

colony analyses information and makes a decision about forager allocation. In the model 

two classes of bee, each with access to different, yet critical, pieces of data interact while 

unloading nectar in such a way that a self-organised value emerges in the form of the search 

time. This value provides both positive and negative feedback, depending on the current 

situation, and so the colony is able to make a wise choice when it decides on its forager 

allocation. The inclusion of receivers in this model provides an inhibition mechanism that 
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is common in decision making models (jSumpterl . 120061 ) and which counters the activation 
provided by forager recruitment. 
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